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Summary 

A nonhydrostatic numerical simulation of a tropical cyclone 
is performed with explicit representation of cumulus on a 
meso-/? scale grid and for a brief period on a meso-y scale 
grid. Individual cumulus plumes are represented by a combi- 
nation of explicit resolution and a 1.5 level closure predicting 
turbulent kinetic energy (TKE). 

The results demonstrate a number of expected and 
unexpected important scale interaction processes. Within 
the central core of the developing cyclone, meso-/i convective 
regions grow and breakdown into propagating inertia- 
gravity waves throughout the lifecycle of the cyclone. In the 
early stages, the amplitude of pressure fluctuations associated 
with the meso-/) scale convection exceed the central pressure 
of the cyclone and strongly modulate its intensity. With each 
meso-/I scale pulsation, the cyclone core increases in strength, 
measured by the central pressure deficit. The increasingly 
strong inertial frequency of the storm core acts to increasingly 
trap the convection induced heating within the core by 
balancing the tangential wind against the low central pressure, 
before the meso-/) scale convection breaks down and sends 
the warmth away as a propagating wave. Eventually, the slow 
manifold’s amplitude exceeds the amplitude of the meso-/? 
scale oscillations and a stable eye region is formed. As inertial 
instability increases, increasingly high thermal warmth can 
be protected in the core, allowing persistent subsidence to 
form and to clear out the cyclone eye. 

On the outside of the eye wall, strong inertial stability 
gradients in the troposphere cause convective warming to 
split the inflow to the eye wall and spawn outwardly 
propagating inertia gravity waves. These waves carry away 
all of the heating forced by convection that is not inertially 
trapped by the eye wall and act as a moderating influence 
on storm intensity. 


Inertia gravity waves are also spawned in the stratosphere 
at the top of the eye wall by the revolution of asymmetric 
cumulus structures. In all instances, the tropospheric waves 
are coupled to the propagating stratospheric waves which 
both move at 35 ms^ 1 , although there are many instances 
where the stratospheric waves seem to have no tropospheric 
counterpart. Hence the anvil top forcing and low level 
breakdown are linked. 

The outwardly propagating inertia gravity waves act to 
initiate outer bands of convection. This initiation is with the 
assistance of low level boundary layer variations of density 
related to previous convection and to virga falling from the 
anvil which moistens and destabilizes the mid levels of 0 e 
minimum. The convection initiated by these waves does not 
move substantially outward with the wave, although may 
appear to develop outward discontinuously. 

1. Introduction 

The notion of CISK (Conditional Instability of 
the Second Kind) was first described by Charney 
and Elliasen (1964) as a way to explain the growth 
of a hurricane from a weak tropical depression. 
The underlying principle is that an ensemble of 
deep precipitating cumulus result in deep tropo- 
spheric heating which lowers the surface pressure 
within the storm, increasing the tangential winds 
of the balanced vortex. The friction layer induces 
radial inflow which, in turn, forces the convergence 
of warm moist marine layer flow into the storm 
to fuel the convection. The critical link in this 
chain of events was the relationship between 
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the surface convergence and the heating function 
which is the role played by deep cumulus convection. 

Early attempts to model this interaction by 
Ooyama (1964) were thwarted by the existence of 
nonlinear instability. Charney and Elliasen (1964) 
were successful in achieving vortex growth because 
they introduced the concept of cumulus param- 
eterization which attempted to represent the 
effects of cumulus without explicitly resolving 
cumulus circulations. Their success together with 
Ooyama’s experiences led to massive observational 
and parameterization development efforts in the 
community to accurately define the heating 
function hoping that tropical cyclone genesis 
could be accurately predicted. Parameterization 
theory spread to other areas of meteorological 
prediction and theoretical analysis where cumulus 
were deemed important, such as general circulation 
theory (Arakawa and Schubert, 1974), wave CISK 
theory (Lindzen, 1974; Raymond, 1975; Raymond, 
1976), mesoscale convective systems (Fritsch and 
Chappel, 1980; Kreitzberg and Perky, 1976) and 
so on. In all applications, the use of cumulus 
parameterization schemes to represent scale 
interaction has resulted in little or no increase in 
the knowledge of the processes they were designed 
to represent. As model resolution increases to 
resolve the meso-/? scale, the underlying cumulus 
parameterization assumption of scale separation 
breaks down. 

The development of three dimensional cumu- 
lus/mesoscale models in the mid 1970s (Klemp 
and Wilhelmson, 1978; Cotton and Tripoli, 1978; 
Clark, 1977) led to massive breakthroughs in the 
1980s in the understanding of cumulus scale 
overturning and gravity wave behavior. At the 
same time, scientists were beginning to recognize 
important aspects of the CISK process beyond 
the effects of Ekman pumping which were highly 
dependent upon the structures of the small circu- 
lations within the cumulus ensemble. For instance, 
gravity waves were proposed as a major component 
of spiral bandedness in tropical cyclones (Kurih- 
ara, 1976; Willoughby, 1977; Willoughby, 1978; 
Willoughby, 1979), the degree to which inertia 
gravity waves were spawned by deep convection 
was shown to strongly modulate the cyclone 
growth rate (Schubert and Hack, 1982), and the 
existence of asymmetric wave-like features in the 
cyclone circulation were shown to significantly 
affect radial momentum transport (Holland, 1983). 


In order to build a scale interaction three- 
dimensional model of a tropical cyclone, all 
thermodynamical and dynamical processes must 
be represented simultaneously. Hydrostatic, axi- 
symmetric, and cumulus parameterized models 
are inadequate for this purpose because of the 
built-in biases and restrictions inherent in their 
design. 

Since Rosenthal (1978) demonstrated that 
Ooyama’s (1964) concern of nonlinear instability 
was somewhat overstated and that explicit con- 
vection models on the tropical cyclone scale 
were possible, there has been renewed interest in 
explicitly representing latent heat release rather 
than parameterizing it. With model resolutions 
allowed under present computing constraints in 
three dimensions, it is still not quite possible to 
realistically resolve the true scale of cumulus 
plumes. So the debate becomes one of whether 
the parameterized approach, with assumed subgrid 
scale cumulus ensembles with clouds of assumed 
simplistic structures, is more realistic or whether 
explicitly predicted convective overturning on the 
scale of mesoscale convective systems rather than 
individual plumes, is more realistic in a model of 
the entire cyclone. 

In a relatively short while, that debate will 
become irrelevant because plausible model 
resolutions will dramatically increase so that 
cumulus can be explicitly resolved realistically. 
In this paper a nonhydrostatic tropical cyclone 
simulation will be presented which has sufficient 
grid resolution to explicitly resolve convective 
overturning on the full meso-/? scale, and for a 
shorter period of simulation, a portion of the 
meso-'y scale near the eye wall. Unresolved con- 
vection was represented by a 1.5 level turbulence 
closure scheme. 

The purpose of this paper is to explore the 
types of internal scale interactions which arise in 
a fully three-dimensional non-hydrostatic model 
of a tropical cyclone when cumulus are explicitly 
represented. As a consequence, it is desirable to 
keep the interaction with the initial state and the 
initial environment as simple as possible so the 
results can be easily understood. It is elected 
to simulate the growth of an axi-symmetric 
modified Rankine vortex perturbation in an initial 
undisturbed marine environment with thermo- 
dynamic structure typical for tropical cyclone 
development and without mean wind. 



An Explicit Three-Dimensional Nonhydrostatic Numerical Simulation of a Tropical Cyclone 


231 


The next section will describe the numerical 
model employed briefly and will be followed by 
a section discussing the experimental design. Next 
the results will be presented for a 56 hour simulation 
period. This will be followed by a section which 
attempts to tie the results together to form a 
coherent picture of the scale analysis process as it 
can be modeled today. Finally, in section 6, 
conclusions concerning the implications of these 
results will be drawn. 

2. Numerical Model 

The numerical model used was the Tripoli (1992) 
nonhydrostatic mesoscale model. The model was 
based on the non-Boussinesq quasi-compressible 
dynamical equations. Model thermodynamics 
were calculated by integrating the enthalpy con- 
serving 6 n along with the total water and several 
ice and liquid water precipitating hydrometeor 
specific humidities. The thermodynamic system 
was closed under the assumption of zero super- 
saturation over liquid which also defined diagnostic 
relationships for vapor and cloud water. For 
these simulations, it was assumed that the entire 
grid lies over a body of water and so the model’s 
terrain following coordinate system was not 
implemented. A spherical horizontal coordinate 
transformation was implemented for this simula- 
tion. 

This numerical experiment utilized no cumulus 
parameterization despite its meso-/? scale hori- 
zontal resolution over the first 50 hours of 
simulation. In this way wave scale interactions 
generated by the convective motions could be 
simulated explicitly, although at times these motions 
could significantly depart from realism due to 
deficiencies of the resolution. Some of the limitations 
of this approach were softened by utilizing a 1.5 
level turbulence closure, where turbulent kinetic 
energy (TKE) was predicted using a modified form 
of the closure developed by Redelsperger and 
Someria (1981, 1982), and subsequently used to 
form the vertical down-gradient mixing terms. 
Some modifications were included to represent 
the effects of saturation, ice phase, and precipitation 
loading on turbulence generation. A more 
complete explanation of the closure is given in 
Appendix A. 

The 1.5 level closure enabled the mixing effects 
of cumulus to be transported by the explicitly 


resolved motions. Upscale transport, such as that 
caused by subgrid scale motions transporting 
subgrid scale fluctuations, and transilient turbu- 
lence (Stull, 1989) cannot be represented by the 
TKE closure used here. Nevertheless, it has been 
found that the simpleTKE closure in conjunction 
with partially resolved cumulus updrafts created 
realistic vertically oriented plumes of high TKE 
which significantly enhance explicitly predicted 
vertical transport locally. Also, horizontal advec- 
tion of the TKE acted to spread the effects of 
cumulus mixing downstream, and in particular, 
around the eye wall. This can be significant since 
the time scale of turbulent dissipation was 
approched by the inertial frequency as the storm 
increases in strength. This will be shown in the 
discussion of results below. 

The surface fluxes of moisture, sensible, and 
radiative heat were calculated assuming an ocean 
of uniform temperature 301 K. The surface layer 
parameterization was based on the Louis (1979) 
surface layer, with surface roughness specified as 
a function of wave height (Delsol et al., 1971). 

Horizontal mixing was based on a deformation 
based 2nd order eddy viscosity closure described 
by Tripoli and Cotton (1982) with a background 
fourth order diffusion as described by Klemp and 
Wilhelmson (1978). The 1.5 level closure was 
deemed inappropriate for the horizontal because 
the issues controlling turbulence were more closely 
tied to numerical as well as physical enstrophy cas- 
cade rather than subgrid scale dry and moist con- 
vective processes which control vertical diffusion. 

Radiative transfer in a cloudy atmosphere was 
predicted with the radiation parameterization 
developed by Chen and Cotton (1983). A modified 
form of the Cotton et al. (1986) (hereafter referred 
to as CEA) explicit microphysics prediction scheme 
was employed. The scheme predicts rain, graupel, 
pristine crystals, and snow crystals. Major modi- 
fications made were to: 

(1) Divide the original pristine ice category into 
a snow and pristine category. The modified 
pristine category was assumed to be composed 
of newly nucleated hexagonal plate crystals of 
uniform mass 1.5x10“ 12 kg. This size was 
determined in separate axisymmetric simulations 
of a tropical cyclone with explicitly predicted 
pristine crystal sizes. The new snow category 
was assumed to follow a Marshall-Palmer 
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distribution with slope and intercept derived 
from an explicitly predicted number con- 
centration per unit mass, similar to the number 
concentration prediction of CEA. The snow 
was assumed to have grown from its nucleation 
size and to be somewhat rimed. 

(2) Aggregated crystals were assumed to be part 
of the snow category. Aggregation of snow 
and pristine crystals was represented by 
appropriate transfers of mass to snow from 
pristine crystals and appropriate adjustments 
of the predicted snow concentration. The 
characteristic snow diameter was limited to 
not exceed 0.33 x 10“ 2 m, which was consistent 
with the Cotton et al. (1986) aggregate distri- 
bution slope. 

The CEA model originally grouped both 
nucleated and new crystals together. Since a 
constant size distribution had been assumed, 
massive nucleation at cold temperatures would 
drastically alter the average crystal size and 
would remove all memory of the growth that 
some of the larger crystals had been through. 
Here, growth processes were assumed to convert 
only a given number of crystals at the specified 
mass to the snow category, which itself was 
logarithmically distributed. Hence, new and 
mature populations of crystals would continue 
to exist where massive nucleation occurred. 
This was especially important to the simula- 
tions of cirrus anvils. 

(3) Graupel was represented with a constant slope 
Marshall-Palmer distribution, with a charac- 
teristic diameter of 2 x 10“ 4 m. This value was 
derived from explicit axi-symmetric simulations 
with predicted graupel mixing ratio and 
concentration where it was found to be charac- 
teristic of the graupel size predicted. 

The numerical infra-structure of the model was 
based on a hybrid leap frog (dynamics) and 
forward (thermodynamics) time integration on 
advection with a semi-implicit time split approach 
to representing the high frequency pseudo-sound 
wave terms. An enstrophy conserving leap frog 
advection scheme was used for the integration of 
wind and pressure while a forward 6th order 
Crowley was used for the integration of the 
thermodynamical quantities. 

A two-way multiply nested Arakawa “C” grid 
system was employed, using a form of the meshing 


technique employed by Tremback et al. (1990) 
and Clark and Farley (1984) modified for a hybrid 
time differencing scheme and to accommodate a 
fourth order smoothing operator in the grid 
center. The grid nesting scheme was also modified 
to allow any nest to move along a specified 
trajectory or to move with the surface pressure 
minimum. 

As shown by Tripoli (1992), the model’s mean 
mass field must be specified externally because of 
the quasi-compressible closure produces a pres- 
sure solution unique only to within a constant. 
This constant was determined by requiring the 
mean surface pressure across the outer domain 
to remain constant by making an adjustment to 
the mean exner function over the entire domain 
at every large grid time step. The adjustment was 
made to the exner function so that local hydrostatic 
balance was unaffected by the change. In addition 
to this adjustment, the outward and inward changes 
of mass flux predicted by the external boundary’s 
radiation condition were forced to balance by 
reducing excessive inflow (or outflow) tendencies 
on a percentage basis. 

3. Experiment Design 

The prescribed initial horizontally homogeneous 
basic thermodynamic structure was taken from 
the temperature and humidity profiles observed 
at Kingston, Jamaica, 36 hours before the passage 
of Hurricane Gilbert. Initial mean winds were 
assumed to be zero. 

The grid locations, spacings, lengths, timesteps 
and implementation period appear in Table 1. 
The outer grid was centered along the expected 
path of the modeled storm, while the medium grid 
was centered over the location of the initial 
perturbation. 

The initial mean state was perturbed by a modi- 
fied Rankine vortex described in Appendix B. The 
vortex was defined by its relative vorticity field 
which was set to decrease from 2 x 10“ 4 near the 
center, to slightly negative outside the center and 
then to increase to zero again. The vorticity field 
is adjusted so that the total area weighted negative 
vorticity exactly balances the positive vorticity at 
a radius of 800 km. Hence no net vorticity was 
added to the flow. The stream function field was 
formed from the relaxation of the vorticity field 
which was in turn used to calculate the wind field. 
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Table 1. Two-Way Nested Grid Setup 


Grid 

1 

2 

3 

4 

Period of activation 

0-56 hr 

0-56 hr 

8-56 hr 

50-56 hr 

# Central latitude 

17.5 N 

16.5 N 

16.5 N 

18.5 N 

# Central longitude 

71.0 W 

68.0 W 

68.0 W 

69.7 W 

# Zonal boxes 

64 

60 

64 

60 

Zonal spacing 

60 km 

20 km 

10km 

3.3 km 

Zonal domain length 

3840 km 

1200 km 

640 km 

198 km 

# Meridional boxes 

64 

60 

64 

60 

Meridional spacing 

60 km 

20 km 

10km 

3.3 km 

Meridional domain length 

3840 km 

1200 km 

640 km 

198 km 

# Vertical boxes 

42 

42 

42 

42 

Vertical spacing 

400-800 m 

400-800 m 

400-800 m 

400-800 m 

Domain height 

26 km 

26 km 

26 km 

26 km 

Large timestep 

120 s 

40 s 

20 s 

6.33 s 

Small timestep 

40 s 

13.3s 

6.67 s 

2.11s 


pressure determined from the integration of the 
outer grid. 

The third mesh was not implemented until 8 
hours of integration in order to conserve computer 
resources during the period when the model was 
beginning to build an Ekman Layer. The fourth 
grid, with a 3.3 km horizontal resolution was 
implemented for a 6 hour period "at the mature 
stage of the simulated tropical cyclone. The limited 
time of its use was both because of the high cost 
of its use (8 hours of Cray-ymp for 1 hour of 
simulation), and because the limited horizontal 
extent of the grid required that the eye wall be 
contained within a box 200 km on a side. 

4. Results 

4.1 Evolution of Structure 



Fig. 1. Initial sounding observed at Kingston, Jamaica at 
0000 UTC 11 September, 1988 


The perturbation produced a maximum tangential 
wind of 1 1 ms “ 1 at the surface at a radial distance 
of 150 km from the vortex center. 

Corresponding virtual potential temperature 
and pressure fields were then calculated using the 
nonlinear balance equation. The outer grid per- 
turbation pressure field was first found using 
Neuman boundary conditions and then normalized 
to the observed sounding at the location of 
Kingston, Jamaica (Fig. 1). The potential tempera- 
ture field was determined by assuming hydrostatic 
balance. The medium mesh was dynamically 
balanced using Dirchlet boundary conditions on 


4.1.1 Initial Rapid Growth Phase - Transition 
to Tropical Storm 

The simulation was initiated at 0000 LST. Because 
the initial modified Rankine vortex was specified 
to be nondivergent, the vortex must have developed 
a divergent wind component by frictional processes 
to initiate convective activity. In this section, the 
modeled transition from a balanced 1 and non- 
divergent vortex, to a developing tropical cyclone 

1 “Balanced” refers to the state where the wind and mass 
fields are linked to each other so that the time tendency for 
horizontal divergence vanishes, as dictated by the nonlinear 
balance equation. 
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Fig. 2. Contours of zonal (u) and meridional (i>) wind components at intervals of 2ms -1 over a region 600 km across centered 
near the vortex core and viewed from the second mesh. Simulation time is (a) 5 hr, (b) 7 hr, (c) 8 hr, (d) 9 hr, (e) 10 hr and f) 1 1 hr 
of simulation 

is described. Admittedly, the modeled transition disturbance and large scale forcing. It was however 

was strongly dependent on the initial structure of instructive to look at the modeled process of 

the Rankine vortex, and such transitions in a real transition in order to gain some insight into the 

developing cyclone were likely to be strongly nature of the explicit interaction between a weakly 

dependent on the nature of the initial perturbing inertially stable vortex and relatively strong, and 
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at times dominant, deep convection processes in 
three dimensions. 

During the first 5 hours of simulation, surface 
friction gradually weakened the initial tangential 
wind in the lowest layer causing it to turn slightly 
inward. Figure 2a depicts the variation of the 
zonal and meridional wind components ( u and v 
respectively) at 5 hours. The intersecting u and v 
contours formed a grid at the center of rotation. 
The u and v grid direction, if defined to be the 
direction of the central v component contour 
pointing toward negative u values, indicated the 
relative degrees of convergent and rotational flow. 
The initial u and v grid direction, with only a 
cyclonic rotational wind, pointed directly north. 
By live hours, a convergence was evident by the 
slight veering of the u and v grid direction. 

By 7 hours, it was found that the surface u and 
v grid in the inner vortex had begun to back, 
indicative of the development of divergence in the 
vortex center. This resulted from a maximizing of 
surface convergence at a radius of 250-300 km 
which was roughly equivalent to the radius of 
peak velocity specified initially for the Rankine 
vortex. This forced a shallow upward circulation 
at the radius of maximum convergence, which 
subsequently forced mass inward, above the surface 
layer toward the vortex center. That, in turn, 
caused the surface pressure in the center to rise 
slightly, reducing the horizontal pressure gradient 
at the center. As a result, the flow inside the 
convergence ring became unbalanced and began 
to diverge outward, strengthening the convergence. 

As the convergence sharpened, a circular ring 
of condensate formed from air being forced upward 
to the lowest condensation level (LCL). Peak 
vertical motions grew very slowly at this time 
before the moisture reached the level of free 
convection (LFC). By 8 hours, Fig. 2c showed an 
amplified convergence ring. As a result of inward 
angular momentum transport, peak winds had 
increased slightly and regained most of the speed 
lost to friction earlier. 

Shortly after 9 hours (Fig. 2d) the forced cloud 
ring reached the LFC and the vertical velocities 
began to dramatically increase from magnitudes 
of centimeters per second to meters per second. 
The deepening convection warmed the air at mid 
to upper levels, and resulted in a more rapid 
lowering of the surface pressure which greatly 
enhanced the surface convergence into the ring 


and hence the vertical growth rate of the ring as 
well. 

In general, because convection favors the 
smallest scales in both the radial and tangential 
directions, any asymmetric perturbation will likely 
initiate a tangential variation of the convection. 
In this otherwise symmetric initial perturbation, 
the only asymmetry, besides the beta effect, was 
the tangential variation of the finite differencing 
truncation error. This had a distinctly tangential 
wave number four variation due to the rectangular 
grid geometry. This very small tangential per- 
turbation grew rapidly and by 10 hours (Fig. 2e), 
a distinct pattern of four dominant convective 
centers appeared, with some weaker meridional 
variation of intensity superimposed from the beta 
effect. 

The role of convection in forcing surface con- 
vergence became dominant over the frictional 
effects by this time. As a result, the surface velocity 
fields converged into these convective systems to 
the point that the flow in the major vortex center 
actually became anticyclonic while the four con- 
vective centers each developed closed circulations 
of their own. In fact, at 10 hours, the surface lows 
beneath the convective systems reached a pertur- 
bation pressure (perturbation from the horizontal 
mean) of — 5.5 hPa while the central pressure of 
the parent vortex pressure was — 3.5 mb. This was 
0.5 hPa greater than the initial surface pressure. 
Hence, the initial modified Rankine vortex had 
broken up into four separate meso-/? scale con- 
vective systems. 

Simulated precipitation formed within the 
individual convective systems and, as with any 
unbalanced convective system, maturity was 
reached and the system progressed into a break- 
down phase (Tripoli and Cotton, 1989). This 
resulted in the formation of deep internal waves 
radiating from the center of the collapse. Phase 
speeds were calculated from animated sequences 
of the vertical motion field and found to be about 
35 ms -1 . The collapse also resulted in the rise of 
the surface pressure beneath the convective centers. 
Convergence, forced at upper levels within the 
parent vortex center, resulted in a strong surface 
pressure decrease to — 1 1 hPa by 1 1 hours (Fig. 2f ). 
Also evident were outward propagating waves in 
the u and v fields indicated by the wave like 
structures in the contours away from the center. 
Three dimensional visualization of these waves 
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Fig. 3. Graphs of (a) peak updraft, (b) peak wind 
speed and (c) minimum surface pressure over the 
storm integration volume of the finest grid activated 
as a function of time from 0 to 56 hr 


(not shown for this particular time) showed these 
waves to be spiral bands of vertical motion, not 
associated with deep convection outside the inner 
convergence ring. The bands sloped upward and 
outward, indicative of upward and outward 
propagation (Durran, 1981). 

It must be recognized that the precise number 
and sizes of the mesoscale convective systems 
formed were probably greatly affected by the 
initiating mechanism from which they were created. 
In this simulation it was the tangential variation 
resulting from numerical truncation that caused the 
four cell convective system structure. The model 
grid resolution of the simulation was also an 
important factor, but not thought to be the 
dominant one since the mes o-fi scale systems were 
over 15 grid lengths across and so not near the 
smallest scale possible. It is most reasonable to 
expect that in a more realistic scenario, the con- 
vection would have developed a strong tangential 
variation as it did here, although the size and 
number of systems would have been dependent on 
variations in initial conditions. It is also likely that 
given no radial variation initially, there is an 
optimal size and number of convective systems 


that would form. Nevertheless, these results 
suggested that even minor departures from axi- 
symmetry would take precedence over any tendency 
for an optimal size in determining the scale of the 
initial convective systems forming in the vortex. It 
is perhaps unsettling to this study that there is no 
reason to believe that this would not have a 
significant impact on the evolution of the cyclone. 

After the 11 hour integration time, another 
period of meso-/? scale growth ensued. Here, the 
low pressure again became most intense beneath 
three primary convective systems, reaching 
— 7.5hPa between the most intense system, 
compared to —6.5 hPa beneath the parent vortex 
central pressure. Note, however, that the parent 
vortex had closed the gap between its intensity 
and that of the meso-/? scale convective systems 
of the first growth period. This was a result of the 
growth of the slow manifold 2 in response to long- 
term inertial adjustment to the convective heating. 


2 The slow manifold here reverses to the balanced part of 
the system, which evolves with inertial frequencies as opposed 
to fast manifold modes which evolve with frequencies related 
to the Brunt Vaisala frequency. 
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The growth and decay cycle of mesoscale con- 
vective system (MCS) activity surrounding the 
cyclone core continued for several hours. During 
this time, the amplitude of the slow manifold 
increased in strength and eventually reached the 
point where it could maintain inflow over and 
above the modulated inflow and outflow of the 
surrrounding convective systems. This occurred 
between 16 and 18 hours of simulation. This 
pulsating development process was also shown by 
Tripoli and Cotton (1989) to be an important 
component of the organization of orogenic 
mesoscale convective systems and could be viewed 
as the realization of the geostrophic adjustment 
process described theoretically by Schubert and 
Hack (1982). 

During this period of growth which began at 9 
hours of simulation, a rapid acceleration of the 
tangential wind occurred. Figure 3 shows the 
evolution of peak vertical motion, peak horizontal 
wind gust, and lowest surface pressure. The peak 
wind gust was somewhat less related to the 
tangential wind during the early stages because of 
the strong divergent wind in the early stages of 
storm development. Sustained tangential winds 
were found to be near 15 ms -1 at 11 hours and 
20 ms -1 at 13 hours and 22ms -1 at 14 hours. 
Hence, the tropical storm phase was reached 
around 14 hours of simulation, approximately 5 
hours after the first free convection. The surface 
low at the center of the storm vortex was — 8.5 hPa 
perturbation pressure of 1002 mb. Since the 14 
hour time was actually in a peak MCS phase, the 
surface pressure was actually 0.5 hPa higher in the 
center than it was at 13 hours. Nevertheless, the 14 
hour time marked the first point where the surface 
pressure remained equal to or lower than the MCS 
pressure lows surrounding the vortex center during 
the MCS phase. 

Figure 4 shows the microphysical, TKE, 9' v , 
vertical motion (viewed from below to see the 
low-level vertical motion), 9 e and potential vorticity 
at 14 hours. The Ertel potential vorticity (Z) used 
is defined as: 

Z = prj-VO, (1) 

where p is air density, rj is the absolute vorticity 
vector and 6 is the potential temperature. The view 
is from the 2nd grid nest in order to attain a storm 
scale perspective. 


Note that at that time, spiral bands were already 
apparent in the cloud water field. Precipitation 
was evident beneath the developing storm. At that 
time, a cirrus anvil was in the initial stages of 
formation. There was a false eye within the con- 
densate field resulting more from the absence of 
convective towers than from forced subsidence. 
This can be seen from the lack of a warm core in 
the 9J. 

The microphysics of the convective regions 
were dominated by pristine crystals at 12-15 km, 
snow at 5-12 km, graupel at 5-7 km AGL, and 
rain below the melting level. Primary precipitation 
growth occurred by the Bergeron-Findeisen pro- 
cess followed by riming within the supercooled 
liquid water to produce graupel. The graupel then 
melted and formed rain. The snow, composed 
primarily of aggregates and rimed crystals reached 
average water contents in excess of 0.6 x 10“ 6 
g m“ 3 at 10 km AGL (above ground level). Graupel 
reached contents near 2 x 10“ 6 gm“ 3 and melted 
entirely before reaching 3 km AGL. The massive 
melting between 3 and 5 km AGL was found to 
have a major impact on the entire storm evolution 
and will be discussed in more detail below. Rain 
reached contents near 5 x 10“ 6 gm“ 3 from colli- 
sion coalescence growth. This produced peak 
rainfall rates of over 9 cm per hour at the surface. 

Figure 4 also shows that the TKE field formed 
vertically erect plumes of vertical mixing in the 
vicinity of the deep convective towers. The TKE 
reached peak magnitudes of 10-20 m 2 s “ 2 , showing 
that considerable vertical mixing assisted the 
explicit vertical transport. Nevertheless, peak 
vertical velocities within explicit convective 
motions were reaching 6-12 ms -1 at that time. 
This was large considering the meso-/? scale of the 
fine mesh. 

Vertically propagating internal waves within 
the stratosphere were evident in the 0' v field. 
Animated visualization of the development showed 
that they emanated from overshooting updrafts 
and were spawned in a spiral with an angular 
frequency at the center equal to the angular 
frequency of the cyclone. This frequency was 
related to the combined lifecycle of the elevated 
cloudtop and the time it took to pass a point while 
circling the storm core. The outward propagation 
phase speed was 35 ms“ \ calculated by clocking 
the wave front through a 40 minute animated 
propagation sequence. The vertical wavelength of 



238 


G. J. Tripoli 



Fig. 4. The simulated three-dimensional storm viewed over the entire second grid (region is 1200 km across) and at 14 hr 
showing: (a) and overhead view of the 0.0 1 x 10 “ 6 g m " 3 condensate surfaces with north pointing up (the rain and graupel surfaces 
are darkely shaded, the cloud water and snow surfaces are lightly shaded, and the pristine crystal surface is partially transparent), 
(b) the view of the condensate surfaces described in a) from the south, (c) the view of the 2 m 2 s 2 TKE surface from the south, (d) 
an oblique view of the 5.1 K 0' v surface from the west, (e) oblique view of the 0.3 ms - 1 (lightly shaded) and —0.3 ms - 1 (darkly 
shaded) vertical motion surface from below with north facing up (note east is to the left here) (f) oblique view of the 0.5 PVU 
surface from the south, and (g) an oblique view of the 355 k d e surface from the south 


the internal waves was about 9 km while the 
horizontal wave length was approximately 150 km. 
Defining the linear gravity phase speed of: 


N 



where N is the Brunt Vaisala frequency of the 
stratosphere calculated to be 2.4 x 10“ 2 and l is 
the vertical wavenumber given by: 



L z 


(2) 


( 3 ) 
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Fig. 5. Same as Fig. 4, except for 26 hr of integration 





and L z is the vertical wavelength, then c p = 34.6 
ms -1 , which is precisely the simulated phase 
speed. 

The vertical motion field, viewed from below, 
showed the geometry of both the vertical motion 
associated with the MCSs and that associated 
with the gravity waves at all vertical levels. Some 
spiral bands were visible aloft in this figure, 
associated with the propagating stratospheric 


waves. Deep vertical plumes of vertical motion in 
excess of 1 ms ~ 1 , not depicted in this figure existed 
in association with the MCS activity. 

The potential vorticity field showed a small 
mid-tropospheric maximum had begun to develop 
in association with diabatic sources within the 
convection at 14 hours. Finally, the 9 e field showed 
growth of the 355 K surface, near the ocean 
surface, but not yet a consistent vertical plume of 
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the 355 K surface on the scale of the 20 km 
averaged grid. The peak 6 e at the surface roughly 
coincided with the strongest surface winds and 
hence the strongest moisture fluxes resulting from 
the bulk surface flux formulation. 

4.1.2 Development of Inertially Stable 
Vortex Center - Transition to Hurricane 

Figure 5 depicts the storm structure at 26 hours as 
the first rapidly deepening phase (which began just 
after 14 hours of simulation) was concluding. The 
local time was 0200 LST. By that time the surface 
pressure had decreased to 950 mb, for a net drop 
of 52 hPa in the 12 hr period since the 14 hour time 
discussed earlier. Peak horizontal winds had 
increased to average speeds approaching 70 ms ~ x , 
clearly reaching the hurricane intensity range. 

The condensate field showed that the convective 
activity had greatly spread outward from the eye 
wall in a banded structure. Low elevation liquid 
only stratiform cloudiness was associated,, with 
much of the banded cloudiness, showing that the 
clouds were forced by the bands rather than the 
clouds forcing the bands. 

The vertical motion field exhibited a good eye 
wall within the troposphere with a series of spiral 
vertical motion bands radiating outward from the 
center at all levels. The spiral bands of vertical 
motion were most apparent within the stratosphere 
and had amplitudes of 0.25 to 0.5 ms -1 . Over 
approximately 30% of the region, however, the 
stratospheric bands extended downward into the 
troposphere, where a second vertical motion 
maximum was found approximately one vertical 
wavelength below the stratospheric maximum. 
The spiral bands were associated with a temperature 
perturbation of 5-10 K amplitude in the strato- 
sphere and less than 2 K in the troposphere. 

Over much of the region, the bands extending 
into the lower troposphere were coincident with 
the condensate discussed above. Animated 
visualization of the condensate field and vertical 
velocity field demonstrated that the condensate 
field was not moving outward with the waves. 
Instead, the waves tended to organize and amplify 
existing cloudiness as they propagated through 
the cloudiness. The clouds within the bands, 
tended to be short lived when compared to the 
propagating bands. This made individual groupings 
of clouds impossible to track for more than a 


couple of hours. Over the period they did exist, 
however, and they seemed to spiral slightly inward, 
with the low level flow. In this sense, they behaved 
similarly to the mesoscale convective systems at 
earlier times. 

The outlying regions of cloudiness were found 
to have a significant impact on the surface layer 
6 e field, creating relatively cold pockets. It was 
found that these cold pockets grew strongest from 
the melting of graupel at the melting level. In fact, 
in at least one case, a 6 e minimum was formed at 
the surface from melting above, which was less 
than the environmental minimum. These theta e 
minima proved to be long-lived and subsequently 
worked their way into the eye wall. Moreover, the 
stronger maximum seemed to initiate new con- 
vection when intersecting with outward propagating 
waves. As these density currents became caught up 
in the vortex, they tended to elongate in the 
horizontally sheared flow and become band-like 
themselves. 

The surface beneath regions of convection within 
the eye wall was found to be associated with a 
relative 9 e minimum itself for the reasons discussed 
above for the outer bands. This produced a mild 
density current and a local pressure increase at the 
surface. This undoubtedly enhanced the con- 
vergence into the eye wall. Moreover these density 
currents together with the density currents merging 
from the outer bands enhanced the asymmetry of 
the eye wall, sometimes resulting in deep convec- 
tion over only half of the eye wall. In fact, it was 
rare to have a completely closed pattern of heavy 
surface rain beneath the eye wall. 

As these density currents moved around the 
eye wall, they seemingly split off the outer edges 
of the surface pressure gradient field, giving rise 
to the outwardly propagating bands which were 
distinctly visible in animated sequences of surface 
pressure. Since the outward propagating bands 
were also closely connected with the stratospheric 
forcing, the split in surface pressure gradient 
occurring at the formation of the band was 
inherently coupled to the stratospheric formation 
associated with the moving elevated cloudtop. 

The predicted TKE structure showed that there 
was a continued existence of deep mixing processes 
in support of the crudely resolved explicit deep 
convection. The TKE within the eye wall had 
begun to form a closed cylinder in the lower 
portion of the eye wall, with individual plume 
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structures remaining aloft. This was indicative of 
the emphasis of storm scale spiraling eye wall 
convection in the lower levels. This resulted from 
increasing inertial stability which suppressed the 
development of divergent wind in response to the 
heating in favor of a rotational flow. Note also the 
extensive anvil level canopy of high TKE produced 
as a result of both shear production from the outflow 
layer and byoyant production from radiative 
cooling at the anvil top. 

The 9' v field showed clear evidence of a warm 
core by 26 hours, forced by subsidence from the 
stratosphere. The development of the warm core 
was associated with the cessation of the growth 
and breakdown process of MCS activity in the eye 
wall discussed earlier. Prior to the development of 
sufficient inertial stability of the vortex, the sub- 
sidence warming forced by individual convective 
plumes resulted in the breakdown of the convection 
in the eye wall due to the destruction of low level 
pressure gradients supporting the MCS inflow. As 
the storm scale vortex increased in inertial stability, 
inertial frequencies increased to the point where 
the tangential wind of the vortex could adjust to 
the warming so that divergent flow from the eye wall 
into the eye would not form and act to destroy eye 
wall convection. Hence by 18 hours of simulation, 
the individual plumes ceased their periodic break- 
down in the eye wall, and the growth of a persistent 
warm core commenced. 

This was shown clearly by Fig. 6 which depicts 
a series of trajectories at the early stage of storm 
development centered around 14 hours and the 
maturing stage centered around 26 hours. Note 
that in the early stage, the updraft trajectories 
spiraled into the center at the surface and then 
moved nearly straight up in deep convection 
plumes. At the later time (Fig. 6b), the trajectories 
spiraled upward in the lowest 5 km and then 
moved vertically upward. As the storm increased 
in intensity toward 50 hours (Fig. 6c) the spiral 
became deeper with little or no explicit vertical 
movement suggestive of vertical convection. Instead 
of spiral convection dominated, which was indi- 
cative of the dominance of the inertial frequency 
over the convective and gravity wave frequencies. 

The continued growth of the inertially stable 
vortex was also evident by the formation of a 
strong potential vorticity maximum of up to 2 
PVU at the storm core. A strongly varying potential 
vorticity field was induced in the lower stratosphere 





Fig. 6. Oblique view of air trajectories through storm core, 
viewed from the south centered at (a) 14 hr, (b) 38 hr and 
(c) 48 hr of integration time 


by the penetration of convection from below. The 
potential vorticity structure at the tropopause 
was characterized by a core of high positive 
potential vorticity surrounded by regions of negative 
potential vorticity induced by convective momentum 
transport from below. 

At this time, the 355 9 e surface had formed a 
vertical conduit within the storm core and from 
the surface to the stratosphere. The increased 
surface values of 9 e were distinctly banded, and 
associated with the banded structure of the 
cloudiness and vertical motion field. 

The increased inertial stabilization of the eye 
wall was clearly evident in the evolution of the 
surface pressure and wind vector fields displayed 
in Figs. 7 and 8. At 14 hours of simulation, 



Fig. 7. Surface pressure field contoured at intervals of 1 hPa every four hours beginning at 14 hr of integration with (a) and 
ending at 54 hr with (k). The pressure at 56 hr (1) is shown 


multiple pressure centers were found, while the 
velocity field was characterized by a series of 
minor vortices within the primary circulation. By 
18 hours a dominant storm core low pressure had 
emerged and by 22 hours a nearly axi-symmetric 


low pressure core was evident. At the same time a 
convergent eye wall flow became more visibly 
consolidated into a circular convergence band and 
by 26 hours it became reduced in size to approxi- 
mately 150 km in diameter. 
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Fig. 8. Same as Fig. 7 except wind vectors are shown instead of pressure at surface level 


Outside the storm core, the tangential variation 
of the pressure structure remained. It was asso- 
ciated with the movement of gravity waves and 
with regions of convection where density currents 
were formed. With each major wave node emana- 
ting from the core, the rate of pressure fall was 


momentarily reduced, suggesting that the gravity 
wave activity was weakening the warm core, and 
hence reducing the rate at which the cyclone could 
deepen. Thus, the convection heating was balanced 
inertially within the core, but continued to form 
wave activity outside the core where inertial 















Fig. 9. Same as Fig. 4, except for 38 hr of integration 


stability was less. This is again consistent with 
adjustment theory posed by Schubert and Hack 
(1982) (Schubert and Hack, 1982). 

4.1.3 Inertially Stable Vortex Attains 
Quasi-Steady Balanced State 
at 38 Hours - Hurricane Stage 

The three dimensional storm structure at 38 hours 
(1400 LST) is depicted in Fig. 9. At this time the 


minimum surface pressure had fallen to 942 mb, 
or another 8 mb since the 26 hour time. In fact, all 
of the pressure fall was before 30 hours when the 
sun rose above the horizon. Overall, the pressure 
evolution was much less dramatic than the previous 
12 hours and, the times of maximum pressure fall, 
roughly coincided with darkness. 

At the 38 hour analysis time the cirrus anvil, 
composed of pristine crystals with contents of 
about 0.5 gm -3 , had expanded beyond the second 
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grid. The vertical motion continued to maintain 
a circular eye wall, however, the banded structure 
outside the eye wall had become somewhat 
fragmented. This seemed to be a result of the 
convection in the outer bands creating gravity 
waves of its own, interfering with the more coherent 
waves emanating from the eye wall. 

Systematic warming of the anvil top in the 
afternoon hours was evident in the 0' v field. It was 
also evidenced by a weaker TKE maximum at the 
upper surface of the anvil and the more shallow 
penetration of the TKE field due to increased 
stability aloft . The lesser vertical penetration of 
convection was also evident in the potential vorti- 
city field, which seemed to produce less excitation 
of stratosphere turbulence than earlier. 

Significant outlying convective activity remained, 
as evidenced by the condensate fields. Comparison 
of this numerical simulation with a host of previous 
similar simulations, suggested that much of the 
outlying convection was supported by the anvil 
itself. It has been found that the amount of outlying 
convection was quite sensitive to the amount of 
water carried out from the center in the anvil. The 
precise reason why has not been investigated at 
this time, however, it is speculated that the precip- 
itation may have moistened and destabilized 
the low to mid levels of the troposphere, releasing 
the thermal “cap” imposed on outlying convection 
created by forced subsidence. This enabled the 
spiral inertia-gravity waves emanating from the 
center to initiate convection more readily. In fact, 
in one simulation similar to the one described in 
this paper except that there was no radiation and 
there was a slightly different setting on pristine 
crystal size, the anvil size was dramatically reduced 
and almost all of the outlying convection was 
eliminated. This resulted in long term dramatically 
coherent stratospheric bands over the entire 50 
hours of simulation. It is likely that result was at 
least partially a result of the radiation difference 
and remains a topic for future research. 

Associated with inertial stabilization at 38 hours 
was the substantial weakening of peak convective 
updrafts to 6- 10 ms “ 1 . This was likely due at least 
partially to the inertial stabilization of the vortex, 
although outlying convective updrafts would not 
be expected to decrease for that reason. Another 
probable reason was the diurnal stabilization of 
the upper troposphere in response to daytime 
warming, which has already been suspected to 


have weakened the TKE and potential vorticity 
fields. The fact that the peak updraft and pressure 
falls increased again after the next sunset, seem to 
support this hypothesis. 

4.1.4 Explosive Deepening Resumes 
After Sunset 

Shortly after sunset, which was at 42 hours of 
simulation (1800 LST), strong deepening of the 
storm resumed and the central pressure fell from 
952 mb at 45 hours to 922 hPa at 53 hours for a 
total fall of 85 hPa from the initial vortex speci- 
fication. 

The simulation at 50 hours was depicted in 
Fig. 10. At this time a massive “tree trunk” of 
355 K 6 e rose upward within the eye wall. Con- 
vection was more limited to the eye wall and inner 
bands. 

TKE field again began to penetrate into the 
lower stratosphere as they did at the 28 hour 
analysis time, although there was an absence of 
strong TKE along the upper surface of the cirrus 
anvil, suggesting increased stability over the 
previous night. It is expected that this may have 
been a result of the accumulated warming of the 
region during the day, the effects of which remained. 

The warm core remained strong, but was con- 
nected to strong warming spread along the base 
of the anvil as a result of longwave absorption of 
upwelling radiation from low levels. The potential 
vorticity maximum continued to be maintained 
also at 50 hours, although greater perturbation 
of the stratospheric potential vorticity field again 
commenced in response to the deepened convection. 

The eye wall diameter remained at 120 km, with 
an eye diameter of about 30 km. Peak tangential 
winds increased in response to the new growth to 
80 ms -1 by 50 hours of integration. 

4.1.5 Extended Simulation with Meso-y 
Scale Resolution 

Because the simulated eye wall had stabilized at 
a mean diameter less than 200 km, it was deemed 
deemed feasible to add a fourth grid with 3.3 km 
horizontal resolution to better resolve eye wall 
processes. The grid was added at 50 hours, and 
the integration was extended to 56 hours. The 
vertical motion field shown in Fig. 10 shows that 
the peak vertical velocities responded immediately 
to the increase resolution but then tended back 
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Fig. 10. Same as Fig. 4, except for 50 hr of integration 


downward. This probably resulted from a short 
term imbalance of the core resulting from less 
overall horizontal dissipation in the balance as a 
result of less numerically related smoothing. The 
peak surface wind speed also momentarily gusted 
to over 85 ms' 1 , but then weakened to previous 
levels. Surprisingly, the minimum surface pressure 
continued on the same trend downward it began 
at 45 hours and then stabilized at between 922 
and 925 mb. 



Figure 11 shows the predicted three-dimensional 
fields after 56 hours. Overall, they did not look 
substantially different than they did with the 
10 km resolution when interpolated to the second 
grid. Some finer structure of the inner bands and 
within the TKE field was apparent. 

To better depict the structure of the inner core, 
selected fields of the simulated eye wall within the 
4th grid were displayed with full resolution. 
Figure 12 depicts the cloud, TKE, w, potential 
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Fig. 11. Same as Fig. 4, except for 6 hr of integration 


vorticity and 6 e fields on the high resolution grid. 
Here the fine structure of the spiraling vertical 
motion field was readily apparent. The eye wall 
was constructed of several asymmetric bands of 
vertical motion. The cloud field was shown to 
almost totally lack cumulus structure and instead 


gave the appearance of a spiral shaped stratiform- 
like cloud. The spiral structure appeared strongly 
in the upper surface of the cirrus canopy, as was 
commonly observed by satellite. 

Note also the more continuous spiraling ap- 
pearance of the 9 e surface maximum. This was a 
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Fig. 12. The simulated three-dimensional storm structure viewed over the fourth grid (region is 198 km across) and at 56 hr 
showing: (a) an overhead view of the 0.01 x 10 " 6 gm“ 3 condensate surfaces with north pointing up (the rain and graupel surfaces 
are darkly shaded, the cloud water and snow surfaces are lightly shaded, and the pristine crystal surface is partially transparent), 
(b) oblique view of the condensate surfaces described in a) from the south showing wave like texture to upper cloud surface, (c) 
oblique view of the 2 m 2 s~ 2 TKE surface from the northwest, (d) overhead view of the 0.5 ms - 1 (lightly shaded) and -0.5 ms“ 1 
(darkly shaded) vertical motion surface from below with north facing up (note east is to the left here) (e) oblique view of the 0.5 
PVU surface from the south, and (f) an oblique view of the 359 K 0 e surface from the south 
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result of the higher resolution depiction of the eye 
wall density current next to the inwardly spiraling 
flow. 

The TKE field was confined to a considerably 
smaller region of the plume. This might be expected 
because a greater portion of the updraft was 
resolved explicitly and so the subgrid scale became 
relatively laminar, except for some local instabi- 
lities within bands. 

The potential vorticity field shows considerably 
more structure presented on this grid. Contour 
cross sections through the field at low levels (not 
shown here) show that the potential vorticity was 
layered into the main center forming a spiral itself. 

The inertial stability can be defined as: 

/MC + /)(f + /) (4) 

where £ is the relative vorticity, / is the coriolis 
parameter, R is the distance to the vortex center, 
and V T is the tangential wind of the vortex. 

Examination of the low level inertial field (not 
shown) showed that the minimum of potential 
vorticity within the spiral was associated with a 
minimum of inertial stability just outside the eye 
wall. In fact, the inertial stability formed a spiral 
relative to the near spiral bands and even became 
negative within the split where the band formation 
was occurring at 56 hours. The inertial stability 
was maximum within the band, but still very 
strong within the eye itself. Hence where outward 
propagating bands were forming, a strong gradient 
of inertial stability was apparent across the region 
of convection. 

This explains why the outer edge of the eye wall 
convection split forming the band. On the inside of 
the eye wall, the subsidence warming by the con- 
vection resulted in gradient wind adjustment 
and thus not weakening the eye wall. On the 
outside edge of the eye wall, angular momentum 
transported in and brought around the storm led 
to weakened inertial stability. This caused 
convergence into the eye wall to be weakened, 
throwing the eye wall convection into a breakdown 
mode. This in turn caused the surface pressure in 
that region of the eye wall to rise while the pressure 
outside the eye wall had fallen. This weakened the 
updraft and gave rise to an outwardly propagating 
inertia-gravity wave. 

It is clear that a coherent relationship between 
the stratospheric wave and the tropospheric wave 


was maintained throughout the inertia-gravity 
wave spawning process. The extent to which 
resonance helped determine the frequency of con- 
vection, the optimal storm rotation rate and a host 
of other possibilities, will remain an intriguing 
topic for future research. 

5. Discussion of Scale Interactions 

The above results described several scale inter- 
action processes occurring during the simulated 
tropical cyclogenesis. These processes have not 
been widely discussed before, primarily because 
previous numerical simulations have used 
simplified hydrostatic, axisymmetric and cumulus 
parameterized frameworks which poorly represent 
such processes. Although all of the processes 
described have not been proven to exist in nature, 
their existence in this simulation raises interesting 
questions and gives valuable insight into how 
three-dimensional tropical cyclones might explicitly 
interact with convection. The following is a brief 
summary of the major types of scale interaction 
noted and their impact on the lifecycle of the 
simulated storm. 

5.1 Simulated Mesoscale Convective Systems 

Convective plumes were simulated to preferentially 
separate into meso-/? scale entities in the early 
stages of storm development, before strong inertial 
balance evolved. These entities, showing both a 
strong gradient balanced and divergent part, 
appeared in the pressure field as Rossby wave-like 
disturbances as described by Challa and Pfeffer 
(1980). The divergent part of the circulation, 
driven by the latent heating aloft, acted to contract 
the vorticity field to produce the wave. 

The simulated embedded MCSs underwent 
growth and breakdown cycles similar to those 
described by Tripoli and Cotton (1989a). These 
cycles were characterized by a period of growth 
where low level convergence was drawn into the 
system at sufficient magnitudes to fuel the system 
growth, while divergence or detrainment was 
produced at upper levels. When low level con- 
vergence was interrupted, the level of no horizontal 
divergence in the updraft moved upward leading 
to the entrainment of lower 6 e air and further 
interruption of the fuel supply. The process has 
been described theoretically by Raymond (1984) 
for the wave-CISK problem, where the growth 
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process was described as an advective mode of 
unstable growth. Raymond found breakdown 
occurred when the surface divergence/convergence 
was upset giving rise to propagating modes of 
inertia-gravity waves. 

The breakdown process, once initiated continues 
because the propagating gravity waves become 
the growth modes. This process was fundamentally 
similar to that which led to the breakdown of 
single cumulus plumes on the 20 minute time scale 
except for the processes which initiate the break- 
down phase. Hence, the field of growing and 
collapsing MCSs continuously emitted inertia 
gravity waves which move about the storm, and 
eventually dissipate or move into the stratosphere 
where the model’s absorbing layer removed them. 

5.2 Gravity Wave Formation Induced 
by Locally Weakened Inertial Stability 

Asymmetric surface deformation fields produced 
by convergence of flow into the developing ring of 
maximum winds produced a tangentially varying 
angular momentum field, which when advected 
around the vortex produced regions of weak 
inertial stability or inertial instability. These 
instabilities locally weakened the convergent flow 
into the eye wall by allowing the flow to 
respond to subsidence warming outside the eye 
wall. This could initiate the local breakdown of 
convection in the eye wall and could give rise to 
the formation of a propagating inertia gravity 
waves in the form of outward propagating bands 
of vertical motion, temperature and pressure. 
These waves, which propagated radially outward 
at approximately 35 ms - 1 , appeared to be always 
coupled with a stratospheric wave of similar phase 
speed. 

The relationship between the upper and lower 
inertia-gravity waves was by no means just a 
coincidence, since the breakdown of the eyewall 
convection at lower levels was also tied to a 
collapsing cloud top. It is interesting that the 
growth and breakdown process seemed to take on 
a continuous form in the eye wall. The convection 
moved around the eye wall with different sections 
of the eye wall containing convections at different 
stages of the process, and thus resulting in the 
continuous emission of spiral gravity wave activity. 
At later times, the simulation seemed to prefer a 
tangential wave number one, with one dominate 


side of upward motion. There were periods however, 
where more continuous eye wall convective activity 
was apparent. 

5.3 Simulated Eye Wall During the Hurricane 
Stage 

As the tangential wind of the vortex increased in 
the hurricane stage, and the core reaches greater 
inertial stability, the tendency for the low level 
wave like MCSs within the eye wall jet were 
suppressed. Whereas individual air trajectories 
move vertically in the early phase of the storm 
within the individual convective systems, the grid 
scale air trajectories in the eye wall move in a 
cyclonic upward spiral during the hurricane stage, 
suggesting a dominance of the inertial frequency 
over the convective frequency within the eye wall. 
Since the predicted TKE magnitudes remained 
high in the eye wall, this motion could be inter- 
preted to mean that individual convective plumes 
may have remained, although perhaps not always 
traversing the entire depth of the troposphere. On 
the other hand, MCSs, or substorm scale groupings 
of convective plumes were eliminated during the 
hurricane stage in favor of a single eye wall system. 

The lowered production of gravity wave energy 
in a hurricane has been discussed before as a 
natural consequence of the geostrophic adjustment 
process (Schubert and Hack, 1982). The lesser 
production of gravity wave energy has been linked 
to greater efficiency of growth and hence an 
acceleration of the growth rate of tropical cyclones 
as they approach maturity. In this simulation, this 
process was seen explicitly as a changing form 
of the mesoscale convective entity from cyclic 
convective entities to a quasi-steady eye wall 
circulation. 

5.4 Simulated Eye During the Hurricane Stage 

Another consequence of the inertial stabilization 
of the eye wall convection was the formation of a 
true “eye”. In the early depression and tropical 
storm stages, an intermittent clearing in the storm 
center resulted by the emphasis of convection 
along the annulus of high winds where the surface 
fluxes of 0 e were highest. This appearance would 
be less likely if a mean wind was considered, which 
would somewhat negate that effect. The early 
eye-like structure exhibited only intermittent 
subsidence within and was too large in diameter. 
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As the storm reached hurricane stage, however, 
and upward motion became more axisymmetric 
along the eye wall, mean subsidence developed to 
the surface producing an eye 20-30 km in diameter. 
This resulted because the inertial frequency of the 
storm increased to exceed that of the formation of 
divergent wind in response to the convective 
heating. As a result, the central warming no longer 
forced nor fueled convective breakdown in the eye 
wall and became long-lived. 

5.5 Microphysics Feedbacks 

It was found that the convective regions produced 
surface 6 e minima, primarily as a result of extensive 
melting in the 3-5 km AGL layer. These surface 6 e 
minima outside the eye wall showed long-term 
persistence and, once formed, tended to be asso- 
ciated with more long term convection. 

The 6 e minima, became entrained into the inward 
spiral of inflow and seemed to be focal points for 
the initiation of new convection by spirally banded 
banded gravity waves passing overhead. 

The ice anvil also appeared to have a substantial 
impact on storm structure and growth. It is expec- 
ted that the radiative feedback was of importance 
because of its affect on increasing upper tropo- 
spheric stability during the day and decreasing it 
during the night. This may have helped to duct 
tropospheric gravity waves beneath the unstable 
anvil top, enabling spiral cloud band formation at 
greater distances from the core. 

The other major impact hypothesized was the 
anvil’s impact on low level stability as a result of 
the virga falling into the 6 e minimum. This would 
have tended to weaken tropospheric stability to 
the point that outwardly propagating inertia- 
gravity waves could more easily initiate clouds or 
unstable convection. 

The present results combined with earlier studies 
such as Kanak (1990), Kanak and Tripoli (1990) 
and other unreported preliminary experiments 
suggest, that in the absence of a high ice water 
content in the cirrus anvil, the diurnal radiative 
effect of the anvil dominates and the existence of 
outlying convective activity took on a diurnal 
fluctuation in intensity. When essentially no anvil 
was present, outlying convection became virtually 
nonexistent. Hence radiative transfer combined 
with microphysical processes seemed to have a 
dramatic influence on storm intensity, structure, 


growth rate and nature of internal scale inter- 
action. 

6. Conclusions 

This paper described nonhydrostatic and three- 
dimensional scale interaction simulations of a 
tropical cyclone. The numerical model successfully 
simulated tropical cyclogenesis using a 1.5 level 
turbulence closure to represent cumulus activity 
on a meso-j 8 scale nested grid system. Although 
the simulation was not of a real storm situation, 
the storm genesis from a prescribed initial vortex 
appeared to be realistic. Important features of 
simulated scale interaction were demonstrated, 
some of which have been implied in the literature 
previously and some of which are apparently new. 
The major conclusions are as follows: 

• The scale interaction within the simulated 
tropical cyclone occurred in two major phases 
summarized as follows: 

(1) The “tropical depression” stage made up the 
first phase of scale interaction, and was 
characterized by the proliferation of meso- 
/? scale growing and decaying convectively 
forced disturbances which spawned deep 
propagating inertia-gravity waves. Waves 
were spawned in both the stratosphere and 
troposphere. The amplitude of the pressure 
ocscillation beneath the MCSs exceeded the 
central pressure of the storm. 

(2) The “tropical storm” stage marked the period 
where the storm’s central pressure exceeded 
the amplitude of the fluctuation in the sur- 
rounding convection. This enabled a 
continuous warm core to grow and maintain 
mean subsidence for increasingly long 
periods. Convective entities within the eye 
wall became longer lived as a result of an 
increasing portion of the induced warming 
being captured within the inertially stable 
eye wall without loss of surface convergence. 

(3) The “hurricane” stage was characterized by 
inertially balanced eye-wall convection which 
maintained continuous eye subsidence which 
allowed the subsidence to work it way to the 
surface. At that time the cyclone pressure 
became nearly axisymmetric inernally 
although important asymmetries existed 
outside the eye wall due to inertia-gravity 
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wave activity. Convection in the eye wall 
became less cellular with updrafts taking on 
the scale and geometry of the eye wall itself. 
The eye wall at times seemed to behave as a 
single convective cell. 

• Simulation with a meso-y resolving grid of 
3.3 km spacing within the eye wall had a surpri- 
singly small impact on the storm evolution, 
although much smoother and more coherent 
looking features were simulated. The results 
suggested that a 10 km horizontal resolution 
was likely sufficient to approximate mature 
tropical cyclone dynamics. 

These simulations crudely begin to unfold the 
explicit scale interaction processes which exist, but 
offer the reader little confidence that the answer has 
been found. Instead, they instill a sense of uneasiness 
that there are important essential processes within 
convective weather systems that have never been 
addressed theoretically or within predictive model 
frameworks. As research numerical models become 
capable of better representing these process by 
simulating a greater range of wave numbers using 
larger grids, a better understanding of the processes 
will be achieved, new questions will be posed and 
requirements for accurate prediction will be made. 
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Appendix A 

Level 1.5 Closure 

The level 1 .5 closure scheme is based on the scheme presented 
by Redelsperger and Sommeria (1981) (hereafter referred to 
as RS) with some modification. For a complete derivation, 
the reader is referred to that paper. Here, the modifications 


to the RS scheme are presented and the final form of the 
implementation of the scheme is given. 

The closure simply involves an alternate procedure for 
computing the mechanical (K M ) and thermal ( K H ) eddy 
mixing coefficients used in the second order diffusion term. 
In the case of the hurricane, the scheme is only applied to 
the calculation of vertical eddy mixing coefficient. The TKE 
is defined as: 


e = 0.5 u' 2 , (Al) 

where u { is the three dimensional velocity tensor and the 
prime refers to the turbulent deviation from the resolved 
scale. Retaining only the largest order terms defined in RS, 
the vertical mechanical mixing coefficient is defined as: 


= C m Le 112 


1 — max(0, C l R i ) 
1 + max(0, C l R i ) 


(A2) 


where l z is the vertical grid spacing, C m ~ -L, and R t is the 
Richardson number, given by: 


N 2 



(A3) 


where N is the Brunt-Vaisala frequency and is the 
deformation tensor. The constant C x is defined to be: 


C,= 


3CC„ 


(A4) 


where C — 4.0 and C 0 =1.2. The thermal eddy mixing 
coefficient is related to the mechanical mixing coefficient by 
the formula: 


KhJC^ 

K m 2C; 

where C u — 4. 

The tendency equation for TKE is given by: 


de 

dt 


Km 


C„L 1 e 2l2 — e 


2 dUi 

3 dx ; 


1 — max(0, CiRi) 
1 + max(0, CiRi) 


(A5) 


(A 6) 


where u ( is a velocity tensor, is the cartesian distance, and 
C £ is the dissipation coefficient which is set to C E = 0.7 except 
near the surface where it is ramped up to C £ — 2.7. 

The buoyancy production (term involving N 2 in A6) is 
the primary term driving the development of deep subgrid 
scale moist convective overturning. The vertical thermal 
buoyancy driven acceleration of velocity can be defined as: 


du 3 \ 
dr , /, 


(A7) 


where 6 0 is a reference state potential temperature and 9 vil 
is a mixed phase virtual potential temperature defined as: 


8-.ll* 61 * 

1.0 + (Ji + q 


(A8) 


where 6 is the potential temperature, q v , q t , and q, are the 
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specific humidities of vapor, liquid and ice in the system. 
This form of virtual potential temperature is useful because 
it takes into account both the effects of vapor and water 
loading in addition to that of temperature in creating buoyancy. 

Following RS, the buoyancy production term is related 
to the vertical heat flux by: 


do 


'= -K g —(^ 
6 0 \dx 3 


L, 

c p T dx 3 




(A9) 


where c p is the specific heat at constant pressure, T is the 
temperature and L m is the latent heat of phase change defined 
as: 


L m = w,L vl + ( 1 + w,)L dU (A 1 0) 

where L vl and L vi are the latent heats of evaporation and 
sublimation respectively and w, is a weighting function used 
to represent the change from liquid to ice in the temperature 
range of 253-273 K given by: 


w t — min { l,max 


(o, 


V 


T- 273.16 
2O0~ 


(All) 


The latent heating term is to take into account the effects of 
heat released or absorbed by phase changes associated with 
either saturated ascent or melting and evaporation on 
unsaturated decent. The vapor change is defined as: 


dq v 


dq m 

E p -mrn{dq m + (q s 


for q v > q s 
q v ),q, + qi) for q„^q’ 


(A 12) 


where q s is the phase weighted saturation mixing ratio, E p 
is an assumed precipitation vaporization efficiency, taken 
here to be E p — 0.5, and q m is the change in vapor mixing 
ratio due to phase change over one grid volume which 
can be estimated (Cotton and Tripoli, 1978) as: 


gAzq v c p T 
RT sL m q v 
Rc p T 2 


(A 13) 


where R is the gas constant and e — 0.622. 

The Brunt- Vaisala frequency implied by Eq. A9, is also 
used to define the Richardson number. This from of the 
Brunt-Vaisala frequency has been found to effectively force 
deep mixing plumes, which represent cumulus reasonably 
well, when the circulation is partially resolved by the grid. 


Appendix B 

Rankine Vortex Initialization 

The initial wind circulation is built from a specified vorticity 
field by solving the Poisson equation for stream function. 
The vorticity field is given by: 

|w,(x 3 )(Ci-C 2 ) for R< R 2 
(0 for R ^ R 2 ’ 

where (, Ci, and £ 2 are relative vorticity values for the 
vortex and two components comprising the vortex 
respectively and R and R 2 = 800km are the great circle 


distances to the vortex center for an arbitrary location and 
the location and the outer limits of the circulation 
respectively. The vortex strength is varied with height by 
the relationship: 


w«(x 3 ) 


j 1 forx 3 <z 1 

j e (-* 3 - Zl / Z2 - Zl ) f orX3>Zl * 


(B2) 


where z^lOkm and z 2 =13km. The inner vorticity 
component is defined: 

Cl=C m axe ,n(01)WJ?l)4 , m 

where R t = 400 km is the radius by which the inner vorticity 
is reduced by one order of magnitude, and ( max = 2.0 x 10 4 
is the center vorticity value. 

The outer vorticity field is added to force the average 
vorticity added by the perturbation to zero and is given by: 


c 2 = sin 



Ki] 



(B4) 


where the horizontal averaging operator is defined for some 
arbitrary variable to be: 


f WA 



where A is the area over which the average is taken, which 
in this case, is the area within R 2 - 
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